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K' 

k 
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Area 

Crack length; half crack length in centrally 
located crack and double edge notched plates. 

Terms of compliance matrix 

Strain -displacement matrix 

Slope of load -displacement curve 

Double edge notched 
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Displacement matrix for element nodes 
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Stress intensity factor 
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Stiffness matrix for the finite element 
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Single edge notched 


A 


u 
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U,V 

W 

r 
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V. 
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iimuN^ 



Distance along J -integral path 
Traction vector 
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Strain energy function for the laminate 
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Load relaxation factor at i^ h ply failure 
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abstract 


The J -integral for quasi-iso tropic (0/+ 45/90)2 S an( ^ 
crossplied (0/90), laminates has been evaluated using 
finite element technique. Constant strain traingular finite 
elements have been used. The program gives displacement 
at all the nodes and stresses and strains in all the elements. 
Stresses are obtained for all the plies and J-integral can 
be obtained for all the laminae as well as for the laminate. 
Three crack configurations namely, single edge 
notched (SECT), double edge notched (DSN)and centrally 
located crack specimens have been analysed,' 

For the formulation of the problem, an equivalent 
stress-strain matrix for the laminate has been defined such 
as to model the laminate by a homogenous and orthotropic 
material, A stepwise method has been used to predict 
l ami nate failure. To account for effect of stacking 
sequence, models with no load relaxation and complete load 
relaxation in failed plies have been investigated. 

The results obtained by present study compare quite well 
with the experimental investigation on SEN specimens. 



CHAPTER I 
INTRODUCTION 

1.1 INTRODUCTION TO FRACTURE OF COMPOSITE MATERIALS 

Composite materials are engineering combination of 
two or more materials to achieve certain physical properties 
not realizable by constituent materials individually. 

Composite materials have high strength and stiffness 
coupled with low density. In addition, they have good 
vibration, fatigue and corrosion resistance, low heat 
conductivity, good electrj cal insulation and favourable 
cost effectiveness. Manufacture of composites requires 
less labour and generates less waste. It has the 
additional advantage of easy p roc es sibil ity. 

Since the composite materials are being increasingly 
used in many sophisticated applications, material 
selection criteria and design procedures are required. 

Global properties like deformation and modulus are well 
understood. Strength criteria, although not well : 

established, are generally agreed upon. So, research 
emphasis is towards understanding of fatigue and fracture 
response of composites. 

Fracture mechanics has now been widely accepted 

as a useful discipline, which deals with inherent macroscopic i 

| 

cracklike defects in materials. It deals with the study j 

of fracture of materials by crack initiation and propagation ! ; 

under various service conditions. It aims at providing | 

| 

guidelines for tolerable size of cracklike defects and j 

[ 

i 
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inspection requirements for avoiding catastrophic failures. 
For metals and ceramics, developments in fracture mechanics 
has led to establishing methods of predicting fracture 
toughness. The fibre reinforced composi f es also seek 
reliable means of characterizing fracture behaviour as a 
part of establishing the fundamental properties. 

Due to heterogeneity of the composite materials 
and complexity of interaction between constituents , 
fracture is accompanied by several processes. Those 
processes include fibre breakage, fibre extraction, 
debonding between fibres and matrix, delamination, fracture 
of matrix etc. In some cases composites fracture during 
stretching along the fibres by breakdown of integrity. 

Even when steady growth of crack is observed, they differ 
substantially from the cracks in traditional materials. 

The problem of characterizing fracture properties of 
composite materials is challenging as the crack propagation 
in them may not be self -similar. D"o to those complexities 
parallel or combined experimental and theoretical analyses 
is required. 

1*2 LITERATURE REVIEW 

In fracture mechanics two basic approaches, namely, 
the energy and stress approaches have emerged. The energy 
approach is due to Griffith [l]. He derived a fracture 
criterion by taking into account total change in energy 
of a cracked body. The stress approach is due to Irwin [ 2 ], 
who relied on stress analysis of a cracked body. 
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Definitive features of fracture in homogeneous 
isotropic solids were identified by Griffith [l] and 
extended by Orowan [ 3 ] and Irwin [ 2 ] for engineering 
materials, Irwin carried out stress analysis by assuming 
the material to be linear elastic. This approach is known 
as linear elastic fracture mechanics (LEFM) . 

Rice [4,5] used a different energy approach 
applicable to linear as well as nonlinear elastic materials 
and under conditions of large scale yielding also. He 
used a two dimensional energy line integral around crack 
tip, now better known as J -integral. This integral 
characterizes the crack tip. Basis of J -integral is 
provided by the works of Hutchinson [6] and Rice and 
Rosengren [7] .. i'fcClintock [3] has demonstrated that 
crack tip stress and strain fields can be described in 
terms of J -integral. The use of J -integral as an elastic- 
plastic fracture criterion has been discussed by Broberg[9], 
Bagley and Landes [10-1 2 ] through extensive experiments 
have shown applicability of J-integral as fracture 
criterion for metals. 

Several researchers applied linear elastic fracture 
mechanics to composite material and found that fracture 
toughness was not independent of initial crack length and 
they used a two parameter LEFM model. Bv a consideration 
of hole size effect. Fuismer and Whitney [13,1+] proposed 
alternative approaches for fracture of laminated composites. 
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These two approaches are point stress criterion and average 
stress criterion. Giare [15] and Agggrwal and Giare [16] 
found notched strength of short fibre composites to be in 
agreement with Whitney-Nuismer criteria. 

The J -integral approach was attempted on glass fabric 
reinforced polyester resins by Smith, Green and Bowyer [17] . 
It has been used to short glass fibre composites by Agarwal, 
Patro and Kumar [18,19]. Srinivasan [20] has used 
J-integral for quasi-isotropic. and crossplied laminates. 

Babu [2l] in a theoretical study has evaluated 
J-integral for short fibre and unidirectional composites 
using finite element method. He used a 2-dimensional 
displacement finite element analysis using constant strain 
traingular elements. J-integral for randomly oriented 
short fibre composites and unidirectional composites was 
evaluated for plane strain and plane stress cases. 

He conducted the analysis for throe different crack 
configuration, namely, single edge notched, double edge 
notched and centrally located crack. His check for path 
independency of J-integral was satisfactory. J for short 
fibre composites are in good agreement with experimental 
works of Patro [22]. 

1,3 3G0PE OF PRESENT WORK 

The present studies arc an attempt to predict the 
fracture behaviour of composite laminates using J-integral 
as a fracture criterion. Symmetric laminates have been 
analysed using equivalent stiffness matrix for the laminates. 
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J -integral has been calculated for individual laminae 
through lamina stiffness matrix. Attempt has been made to 
predict first ply failure and analysis after first ply 
failure has been carried out by modified stiffness matrix. 

The analysis has been carried out for quasi -isotropic 
[0/+ 45 / 90 ] 2s and crossplied [ 0 / 90 ]^ s laminates for three 
different crack configurations, namely, single edge 
notched, double edge notched and centre located cracked.' 

Formulation of the problem is given in chapter 2. 
Finite element solution technique adopted for the present 
study is also discussed in chapter 2. Chapter 3 deals 
with the numerical result s obtai nod. The various models 
adopted to predict fracture behaviour of laminates are 
also discussed in this chapter, in the light of experimental 
findings. Conclusions ana some suggestion for future work 
are included in chapter 4* The computer program used for 
the present study has been given as an appendix." 



CHAPTER II 


J-INTEGRAL AS FRACTURE CRITERION FOR COMPOSITE MATERIALS 


2. 1 INTRODUCTION 


Most of the research work on fracture of composite 
materials, so far uses linear elastic fracture mechanics 
(LEFM) concepts and requires accurate stress analysis at 
the crack tip. Due to stress singularity at crack tip, 
this is difficult even for homogenous isotropic materials* 
In heterogenous composites, the situation becomes more 
difficult due to complex crack tip geometry. 

The J-integral is a parameter which can be evaluated, 
due to its path independency, without focussing attention 
at the crack tip, but can still characterize the crack tip. 
J-integral does not require a prescribed crack extension 
trajectory and is hence more adaptable to characterize 
cracks in composite material. 

The J-integral, as proposed by Rice [4,5] is a two 
dimensional energy line integral defined as 

u i 

(Udy -T^ - ds ) (2.1 ) 

where p is any contour traversing in counter clockwise 

sense and enclosing the crack tip as shown in Fig. 2.1^ 

U the strain energy density is defined as 

e 

t mn 


U = U ( e ) 
v mn 


a ±d d ®iD 


where c . . and e . . are stress and strain tensors 
13 13 
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respectively* T. is the traction vector and is given "by 
1- = a - n . 

■where n^ is the normal vector to p . is the displacement 
vector and s the arc length along p. 

2,2 FORMULATION OF PROBLEM 

For evaluation of -I -integral, stresses, strains and 
derivatives of displacement are to b- found out. To achieve 
this end total potential energy stored in the laminate has 
been minimized by Raleigh -Ritz method [2?]. The finite 
element program, then selves the pro clem with prescribed 
boundary conditions to yield displacement at all the nodes, 
their derivatives, strains and stresses. 

2.2.1 Expression for Potential Energy 

For symmetric laminates subjected to inplanc loads 
only, displacement and strain fields are same throughout 
the thickness of the laminate, end thus donot vary from 
lamina to lamina [24-]. 

Due to different fibre dire ction , stress-strain 
matrix vary from lamina to lamina and honce the stresses 
in them are not same even with same strain field, 

JL1~ 

Row if the stress-strain relationship of i 
lamina in a laminate be [D]j_ and its thickness be t^, 
strain energy stored in it is given by 

U i= 2 X '-'JiMi dT= L-eJi IPH ^(2.2) 

i i 1 
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where L£i is the row matrix Lf-i e 2 ^12 J an & ( e j 
is the column matrix. So the total strain energy stored in 
the laminate will he 

U = z u =| L. e J (£,[ dv ±) \ £ ] (2.3) 

i i 


Now if the laminate is a flat plate 

dv^ = dA t ± (2.4-) 

and eqn, 2.3 can he written as 


U = 




(2.5) 


Now an equivalent stress strain matrix for the 
laminate is defined such that 


[D ] 0 t 


c 


£ [dL q 

i „ 


where [d] is the equivalent stress-strain matrix and 
c 

is the total laminate thickness. 

Equation 2.5 can he rewritten as 


( 2 . 6 ) 


u = J A t_ e J [D] 0 i e i t c <“• 


(2.7) 


Now in the absence of hodv forces expression for the potential 
energy of the laminate can he written as 


I = 


L S J L D 1, 



t 

c 



i u as 
3 3 


( 2 . 8 ) 
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2*2,2 Bonn dary C on di ti on 

Three different crack configurations namely single 
edge notched, double edge notched, and centrally located 
crack plate have been analysed. These plates have been 
subject.:! to mode I loading condition as shown in Fig. 2.2. 
Boundary condition in all of them are 


Along AB 


V = 0, r =0 
’ xy 


Along BC and DA a = 0 ; x - 0 (Traction free) 

x xy 

Along CD V - & (applied displacement) 

and x =0 
xy 

crack surfaces are stress free. 

2.3 FINITE ELEMENT SOLUTION TECHNIQUE 

Expression for the potential energy of the laminate, 
Eq, 2.8, contains only first derivative of displacements. 
So compatibility of u, v and completeness of 


u, v, 


Qjl. 

ax 


d v 
ax » 


d u 

ay ’ 


are required [23]. 


ay 


Traingular elements with three nodes with following 

interpolation function satisfy the above requirements. 

e ( U 1 
{u] = a +■ bx + oy = Ul !J 2 n 3 J j u 2 


Thus 


u- 


* v 0 

jv J = a* + b'x + c'y = [__Nj 


N ? Hi j 3 Vo \ 


ri) 

2 

V- 


(2.9) 


where j__N^ Np N^J are Interpolation functions. 

Substituting equation 2. 9 in 2. 8 and applying the Raleigh- 
Ritz method [23], one gets 



oo 



THREE CRACK 
WITH BOUNDAR 

I 





1 2 


[kl 


a ! 


>(ne) 


= |p] 


(ne) 


( 2 . 10 ) 


e r p T 

where [k] = JJ [b] [d] [b] t di = stiffness matrix. 


A(e) 


[B], 


t>N. 
2 x 

0 


£ 5 l 

S 3 7 


0 


BN. 


c- y 


;>n ± j 

X i 


,(ne) 


and j d\ 


! “"i / 

L\ 


i is the node numbs r. 

So for traingular elements 


[B] = 


bn. 


ax 

o 


BN. 


X 


o> y 


ne 


and i d 


a n„ 


bn 


1 


a y 


bn, 


Bx 


0 


bn„ 


dx 


^ U 1 
V 1 

/ U 2 


u. 






0 


a n. 


a y 


BN, 


d x 


0 


bn 0 bn 


By 3x By 


0 


BN 


a y 


aa 3 
3 x 


( 2 . 11 ) 


( 2 . 12 ). 
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( 

1 


P 




P 


M 


ds = Load matrix 


(2.13) 


Ng = interpolation function along the boundaries of the 
finite element. 

How the specimen to be analysed is divided into a 
finite number of traingular elements. 1 moderately fine 
element grid is maintained near crack tip. The size of 
elements is increased as it is moved away from crack tip 
because of less variation in stress. Element stiffness 
and load matrices are obtained and assembled. After 
applying appropriate boundary condition t displacement, 
derivatives of displacement and strains are obtained. After 
this using the lamina properties stresses in different 
laminae and consequently J is determined for all laminae 
seperately. In view of the path independency, contour 
is chosen along the boundary of the specimen as it is 
convenient. 


2. If A GO ERECT IOH FACTOR 

For homogenous and orthotropic material , relation 
between I and stress intensity factor K is given by [25] 


T 



(2. Ilf) 


t' 


where 



a ll a 22 )l/2 j- ( 
2 


a 22 

a ll 


^ 1/2 + 


a 66 + 2a 12 ] 1/2 


2a. 


11 


(2.15) 



Ilf 


a. . are terms of the compliance matrix of the material. 

-L J 


how if the elastic modulus in loading direction he 


E 


E 


then J = ( — y ) 


'1 


E 


E, 


how from the finite clement formulation only 

E 


factor is evaluated, so a correction factor 
introduced to get J value for the laminate. 


E 


■1 . 
-y- IS 


( 2 . 16 ) 

K 2 


*1 



CHAPTER III 


RESULTS AND DISCUSSION 

This chapter deals with the numerical results obtained 
for quasi-isotropic and crosspiied laminates and their 
comparison with experimentally obtained values. An attempt 
has been made to predict laminate failure by using properties 
of constituent laminae. A stepwise method has been used to 
predict laminate failure. 

3.1 PROCEDURE FOR LAMINATE ANALYSIS 

For a given specimen, J-integral was evaluated at 
different displacements for the laminate as well as for the 
laminae j constituting the laminate. At smaL 1 displacement 
none of the plies fails. Ac the load increases, plies start 
failing. A oly was assumed to fail when the applied 
displacement exceeded a critical displacement for the ply. 

The final laminate fracture is assumed at the critical 
displacement for the laminate obtained experimentally. Tie 
basic lamina properties and the critical displacements for the 
purpose of present investigation were obtained experimentally 
and are given in Tabic 3 » 1 . 

After a ply fails, its contribution to the overall 
laminate behaviour changes. This change is incorporated in 
the analysis through modification of Its stiffness matrix. 
There are different approaches which can be used to modify 
the stiffness matrix [ 24 ]. One approach suggests that the 
stiffness matrix of a failed lamina should be made zero. This 
would mean that the lamina would carry no additional load in 
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any direction. This will be a conservative approach since, 
for example, when a ply fails by transverse tensile failure, 
it may not carry any additional load in transverse direction 
but can certainly carry additional load in the fiber direction 
A more practical approach under such circumstances would be 
to assume that while the transv-rsc tonsils modulus and shear 
modulus of a failed ply arc zero but the longitudinal modulus 
remains unaffected, Tim. latter approach has been used in the 
present analysis. 


TABLE 3,1 

LAMINA ELASTIC CONSTANTS 


Longitudinal modulus, E L 


31.9 GPa 

Transverse- modulus, 

— 

9.7 GPa 

Major poisson’s ratio, V 

= 

0.2^ 

Shear modulus, 


3.2 GPa 

Lamina.tc s 

Critical displacement 
(mm) 

0° unidirectional 


0,88 

90 ° unidirectional 


0,08 

4-5° unidirectional 


0.10 

Quasi -isotropic, (0/+ 45/9° ) 2s 


0 . 75 * 

Crossplied, ( 0 / 90 )^ g 


0.375* 


* From reference 20 
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The J-integral evaluated through a stepwise procedure 
may he represented hy a piecewise continuous curve as shown 
in Fig, 3.1, ’where each discontinuity point represents a 
ply failure. In the segment between two discontinuities the 
J-integral, (J)^, can be written as follows 

(J) ± = + ( * 1) ± (3.D 

Where (J)^_^ is the value of J at the (i-1)^^ 1 
discontinuity and ( A J)^ is the increamental J in i® 1 segment. 
When all the plies are intact, the J-integral is evaluated 
using the original stiffness matrix for the laminate. The 
increamental J is calculated using modified stiffness 
matrix. 

Another point which must be considered at this point 
is what happens to the load carried by the ply when it fails 
That is, when a ply fails what is bhe extent to which the 
load carried by it drops. In a real situation, it may be 
between the two extreme situations, namely, one in which the 
load in felled ply drops to zero ( complete load relaxation) 
and two, the load doesnot drop at all ( no load relaxation). 

It can be easily shown that the later condition would 
show greater increament of J compared to the former condition 
for the same displacement increament. In the present 
investigations, results have been obtained for both the 
conditions, 

3, 2 J-INTEGRAL FOR LAMINATES 

J-integral ms evaluated as a function of applied 
displacement (6) for quasi-isotropic and crossplied 
laminates. The ratio of crack length to specimen width was 



-INTEGRAL 



FIG. 3.1 J vs 6 PLOT FOR A LAMINATE 


19 


varried from 0,1 to 0,7. The following orach configurations 
were used: 

i) Single edge notched (SEN) specimens 

ii) Double edge notched (DEN) specimens 

iii) Centrally located notch specimen 
J ys ^ curves for SEN specimens are shown in Figs. 3.2 and 
3.3 for quasi-iso tropic and crossplied laminates respectively. 

1 is plotted against ^ 2 in Figs. 3.4- and 3*5. Since <1 vs 6 2 
plots are straight lines passing through the origin, it is clear 
that <1 -integral is proportional to the square of -applied 
displacement for all these cases. In fact, this behaviour is 
expected. In the present analysis, it has been assumed that 
the laminate stress-strain ( or load -displacement) behaviour is 
linear. The load -displacement curves for a laminate with 
different crack lengths may be represented as shown in 
Fig, 3*6. For the linear relationship, the potential energy 
stored is given by 

U = | c 62 (3.2) 

where l c T is the slope of P- 6 cruve. It should be 
noted that c is a function of crack length and type of 
crack, but is independent of applied displacement 1 s 1 . 
Therefore J-integral at applied displacement 1 6 1 will be 

r=-|2- =(-1 4^ > 62 = k 62 0-3) 

d. $ a 

where k = - ~ is again a function of crack configuration 

only. 



, KJ/m 2 



DISPLACEMEN 


FIG. 3.2 J vs 6 PLOT FOR 


KJ/m 2 



DISPLACEMENT , mm [ 

FIG. 3.3 J vs 6 PLOT FOR DIFFERENT (a/w) | 



KJ/m 





, KJ/m 



(DISPLACEMENT) 2 , mm 2 
FIG. 3.5 J vs 6 2 PLOT 



2 


I - Crack Length 'a{ 
II - Crack Length 'a 2 
a 2 ^>a^ 



DISPLACEMENT 

FIG. 3.6 LOAD DISPLACEMENT CURVES FOR 
DIFFERENT CRACK LENGTHS. 
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In view of the above discussion, the J- ^ relationship 
can be studied, without any loss of information, through the 
parameter (.3/ 6 w ) representing slope of J vs ^ li nes in 
F. : rs, 3,4 and 3.5. 

For a particular laminate with given crack 
configuration, the family of J -6 curves can be reduced to 
a single curve of (l/ 6^) plotted against crack length to 
specimen width ratio* 

For the purpose of analysis, j/ for the laminate 
was obtained with its primary stiffness matrix. After a ply 
failure, (J - / values were again obtained with the 
modified stiffness matrix. T. 'is process was continued till 
the laminate failed. 

o 

Values of (j/ <5 ) wore obtained for quasi -iso tropic 
and crossplied laminates having three different crack 
configurations, namely (i) SEN , (ii) DEN and (iii) centrally 
located crack, j/6^ values are given in Tables 3.2 to 3.7 
and shown graphically in Fi-~s. 3.7 to 3.-12. T-e values 
have been shown for the laminate when all the plies are 
intact and also for the laminates after initial failure, i.e. 
when plies have failed in ascending order of their critidal 
displacement. These curves are used to predict the overall 
laminate behaviour. It may be noticed ( Figs, 3*2 and 3*3) 
that for SEN specimens J" vs. 6 curves for a/w ratios of 
0.5, 0.6 and 0,7 are quite close to each other for quasi- 
isotropic as well as crossplied laminates. This is in good 
agreement with experimental observations that T-integral curve 
is independent of crack length when a/w >0.5 [18,20,22 1. 
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TABLE -3 , 2 

Single Edge Notched 

Quasi-iso tropic Laminate 

a > J/S 2 in kJ/m 2 -mm 2 

,w . . I 

, ALL PLIES | 90 PLIES , 90 and + 45 ° 

f ' I FAILED PLIES FAILED 

- ** 4 -“ — - 


0.1 

: 9.78 

| 8.48 

6.10 ; 

0.2 

; 32.15 

27.82 

19.65 | 

0.3 

49* 88 

' 43. 27 

29.55 i 

0.4 

1 73. 25 

| 

63.62 

46.22 | 

o.5 

87*55 

: 76.22 

56.32 

0, 6 i 

99.98 

t 

! 87 . 20 

65.42 

0.7 

; 102.0 

89.02 

67.48 

< 



TABLE 3.3 




Single Edge Notched 




Crossplied Laminate 


i 

a I 


j/ S 2 in kj/rn 2 -mm 2 | 

i 

w 

! 

- ' 

ALL : 

PLIES INTACT 

90 PLIES FAILED' 

0.1 


21. 95 

18.' 2 

0. 2 


81.25 

: 71.5 

| l 

| 0,3 


113.35 

94. 35 

[ 0.4 


154. 60 

123.30 

o.5 


170.3 

131.6 

* 

0.6 

I 

181. 25 

! 136.15 

j 

0.7 

f 

[ 

178.55 

1 133.2 

! 



TABLE _3,4 
Double Edge Notched 

luasi-isotropic Laminate 
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2a 


j/ 6 2 in 


w 

ALL PLUS 

INTACT 

90° PLIES 
FAILED 

90° and + 45° | 
PLIES FAILED j 

0.1 

6.83 

5.92 

4.28 

j 

0. 2 

14.27 

12. 28 

8.67 ! 

0.3 

20.88 

17.95 

12.77 | 

j 0.4 

27.18 

23.41 

16.68 

1 0.5 

I 

33.09 I 

28.42 

S 

j 20. 27 

i 0.6 
* 

39.35 ! 

34*06 

! 24* 28 j 

' 0.7 

| 

i. 

46.71 

40.52 

j 29. 00 1 

1 i 

1 — . — a 


TABLE 3.5 

Double Edge Notched 
Crossplied Laminate 


2a | J/ 6 2 in kJ/m 2, - mm 2 


; W j 

1 

1 

ALL PLIES 

INTACT 

. - 

90° PLIES ! 

FAILED | 

. - - - - - j 

| 

0.1 

15.35 

12.51 | 

0.2 

I 

36.27 

32.24 

i 0.3 

51.14 

44.70 

| 

1 0.4 

66.74 

| 58.61 i 

1 | 

! 0,5 

79.53 

) 

i 68. 91 j 

| j 

; 0. 6 

j 94.27 

i 80.96 ! 

i i 

| 0.7 

! 

$ - 

108. 25 

1 91.04 | 

! 1 

t \ 

1 J 



TABLE _3.6 

Centrally Located Crack 
Quasi -isotropic Laminate 
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2 a : 

W ; 

J / in 

kl/m^-mm^ 

j 

1 

\ 

j 


ALL PLIES j 

INTACT j 

90° PLIES 
FAILED 

90 ° and + 45° ! 
PLIES FAILED J 

0.1 ! 

i 

5.90 ! 

1 

5* 12 i 

| 

j 

3.66 | 

C\J 

* 

0 

12. 31 ‘ 

10,66 

7.47 | 

0.3 ; 

18. 51 

16.02 

11.27 ! 

0.4 

; 25.51 

| 

22.08 

15. 56 ; 

0.5 

32.98 

28 . 54 

20. 25 ; 

| 

0.6 

42.59 

1 36.82 1 

26. 34 ; 

‘ ! 

0.7 

54.21 

! 

1 

; 46. 88 

1 1 

33.91 ; 

i 

1 
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TABLE _3* 7 


Centrally Located Crack 
Crossplied Laminate 


2a 

w 


j/ 6^ in kJ/m^-mm^ 


0.1 


ALL PLIES 
INTACT 


14.05 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


33.06 
48.01 
65.29 
80.48 
99.22 
117. 88 




? 

| 

i 


90° PLIES 
FAILED 

12.08 

30.28 

42.94 

58.02 

69.91 

83.97 

96.34 


"~1 



i 

t 

1 

{ 

1 


j 

j 

\ 


) 



(J/6 2 ) , KJ/m 



a / w 


FIG. 3.7 J/6 2 


FOR QUASI -ISOTROPIC LAMINATE 



(J/6 2 ) ; KJ/m 2 


£ 50 
E 


QUASI -ISOTROPIC 
CENTRALLY LOCATED CRACK 
'I - All Plies Intact 
II - 90° Plies Fdfled 
III - 90° and - 45° Plies Failed 



2a /w 

FIG. 3.9 J/6 2 FOR QUASI -ISOTROPIC LAMINATE 


(J/6 2 ) ; KJ/m 





(J/6 2 ) ; KJ/m 2 


140p 

CR0S5PLIED (DEN) 

120- 1 ‘ A| l Plies Intact 

, II “ 90° Plies Failed 

si 

| 100 - 



2a/w 


FIG. 3.11 J/6 2 FOR CROSSPLIED LAMINATE. 


(J/6 2 ) ; KJ/m 



2a/w 


FIG. 3.12 J/6 2 FOR CROSSPUED LAMINATE 
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Ehe procedure adopted for evaluating J-integral of 
the laminate through the stepwise procedure can he explained 
as follows- . Figure 3.13 shows an idealized load displacement 
curve of a specimen. The path OABC * represents the model in 
which no stress relaxation takes place when a ply fails. 

The? path OADEFG represents the model assuming complete load 
relaxation upon a ply failure. For the analysis hero, a 
more practical case where a partial load relaxation takes 
place ( represented Toy OAHIJK ) has been assumed. The 
slopes c-^, C 2 and c^ of the load-displacement curve segments 
between ply failures depend only on the initial crack 
length and are independent of the extent of load relaxation. 

If tk denotes a load relaxation factor at a ply failure it 

can be defined as follows ( with reference to Fig. 3*13) . 

- JI 

n l AD 

^2 = ~ir <3.1**) 

To obtain J-integral at a given displacement 1 & l , area 
under the load displacement curve is first obtained . The 
throe areas shown shaded in Fig. 3 .13 can be easily 
obtained as 

A i = k e i 

^2 = I °2 (62 - 6 1 5 

A 3 = (l-T^) (^-Cg) S p( M -| ) ^3r5) 

So the expression for potential energy can be written as 
u = a x + A£ + a 3 = \ c x &*+j (6 2 -6p + (i-n 1 )(c 1 -c 2 )6 1 

( 3 . 6 ) 




6] 6 62 

DISPLACEMENT 


FIG. 3.13 LOAD -DISPLACEMENT CURVE FOR A LAMIN^ 
WITH PLY FAILURE. I 
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Now 'J-inborr.iL is given by J 


a u 

a ? so ^Liat 


J 


1 


d c 
da 


Sc, 


= 1 - I -a-V 


h O , X C_ 

(--i _ k_2 ) 6, ( 6.5,) 


a 


d a 

1 ^ 


(3.7) 


Defining ^ | and k„ = - i -—2 


da 


2 2 da 

2\ ■ , 


whore k^ and k^ are. the values of (J"/ 3 ) in those segments, 
J= + k 2 (d 2 _6^) +2 (1 _n i )(k 1 _k 2 )d 1 (G _ 4 1 ) 

= k 2 6 2 + (k 1 -k 2 )d 2 +2(1 -T) 1 )(fc 1 -k 2 )S 1 (§. 

x* 4-Vj 

So if 1 0 1 is in n segment then 
n 


(3*8) 


J=V 2 +Z [ Oc^, -k ± ) 4 ( , + 2(l J l 1 _ 1 )(k 1 . 1 - 

(?.9) 

Equation 3-9 has been used to obtain J- 6 plot for the 
laminates. 


3.3 COMPARISON WETTI THE. EXPERIMENTAL RESULT ! 

I 

Predictions of J c for SEN specimens could be compared I 
with the experimental results which are available only for ; 

j 

SEN specimens [20], Sr inivasan [20] found that J vs 6 
curves for quasi -isotropic as well as crossplied laminates t 
arc independent of crack length when ^ ^0,5. The 
experimental value of J c for a laminate is obtained from this 
common J- 6 curve » In the present study, as has already been 
pointed out, <J- 6 curves for SEN specimens, though change 
with ovack length, aro quite close to each other for 
a/w > 0,5. 
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In .]'< ir ’(\senfc study .1 values for the laminates 
has been obtained for three different models namely (i) 
when all tV plies fail simultaneously (ii) when plies fail 
successively with no load relaxation in failed plies , and 
(iii) when plies fail successively with complete load 
relaxation in failed plies* 

J vs curves for these throe models in the range 
0.5 <c a/w ,< 0.7 arc shown in Figs. 3*14 to 3.19 for quasi- 
iso tropic and crossplicd laminates. Values of J obtained 
through different models are given in Tables 3,8 and 3*9 
alongwith the experimental values. It is observed that the 
experimental and theoretical curves are quite comparable. 
However, the model in which all plies fail simultaneously 
predicts a higher value of J than the experimental value. 
This is expected because the off-axis plies do fail before 
the final fracture of the laminate. This model doesnot 
require a step-wise prediction method to be applied and 
hence is simple to apply. 

The models in which plies fail successively predict 

lower values of J" c than the experimental. But the model 

which assumes no load relaxation in failed plies predicts 

j- & curve and J value closes to the experimental ones. 

c 

The difference between the predicted values of J c through 
the models with no load relaxation and with complete load 
relaxation is about 10$ and as may be expected former model 
predicts a higher value of J Q , Thus it may be concluded 
that it is basically the extensional s biff ness matrix of the 
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QUASI- ISOTROPIC (SEN) 
Successive Ply Failure 
No Load Relaxation 
° Experimental 

I a/w = 0.5 

II a/w = 0.6 

III a/w =0.7 
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DISPLACEMENT, mm 

FIG. 3.15 J vs 6 PLOT FOR QUASI - ISOTROPIC LAMINAT 
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FIG. 3.16 J vs 6 PLOTS FOR QUASI - ISOTROPIC LAMINATE 
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FIG. 3.18 J vs 6 PLOT FOR CROSSPLIED LAMINATE. 
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FIG. 3.19 J vs 6 PLOTS FOR CROSSPLIED LAMINATE 
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0,5 

0,6 

0.7 


TABLE -3. 8 

Quasi -isotropic Liminate 


J c in 

Present 

kJ/m 2 

■work 

- 

j "ExperfL" 

| mental 

Simultaneous 
ply failure ; 

Successive ply failure 

no load Complete load 

relaxation' relaxation 

! 

I » 

i [20] 

1 

1 

1 

^ 9. 25 

! 

35.75 | 

31.95 

f 1 

i 

! 

i 

56, 24 

4 1. 30 ; 

37.10 

46.5 

57.37 

4 2.46 ; 

t 

38. 25 | 



TABLE _3*9 
Crossplied Laminate 



J 0 in 

kJ/m^ 




Present 

work 


Expenm 

mental 

a/w 

Simultaneous ! Successive 

1. ...... . . 

ply failure 

1 1 

O 

eg 

1 — 1 


ply failure 

no load 

complete load 

! 

( 

Relaxation 

! 

I 

1 

relaxation 

i 

s 

i 

I o.5 

i 

- - j 

23.95 

20.58 

18.76 


| 0.6 

25.49 

21.56 

19.44 | 

22.65 j 

[ 

O 

• 

-o 

25.11 

21. 16 

19.03 

' 

f 


* 
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laminate that controls its fracture toughness and stacking 
sequence may not influence the fracture toughness signifi- 
cantly. The lower value of J predicted by step-wise 
procedure may be attributed to the fact that in these models, 
- T of the failed plies have been assumed to be 

zero and so in effect those layers carry no additional load 
in transverse direction. However, in a practical situation 
these layers do carry additional transverse load away from 
the dominant crack. 

It may be pointed out that the final fracture of the 
laminate was assumed to occur at the experimentally obtained 
critical displacement of the laminate and not at the critical 
displacement of the 0° plies. Although this limits the 
applicability of the step-wise procedure , it is a necessary 
assumption in view of the experimental observations. 

Actually this is a consequence of a change in the critical 
displacement of 0° plius due to prior failure of off-axis 
plies and the resulting stress concentrations. 



CHAPTER IV 


• CONCLUSIONS AND SUGT-ESTION FOR FUTURE WORK 

I-integral for quasi -iso tropic (0/+ 4-5/90 ) 2 and 
crossplied (0/90)^ g laminates has been determined by finite 
element technique. 

A stepwise analysis procedure was employed to obtain 
Jl. After a ply failure stiffness matrix of the failed ply 
was modified for further analysis. Models with no load 
relaxation and complete load relaxation in failed plies 
were investigated. 

Based on the results presented in Chapter III the 
following conclusions are made. 

(i) U~ <5 curve obtained for 0.5 .< a/w ^ 0,7 for SEN 
specimens are quite close to each other. This compares 
favourably with experimental finding of an unique 

J- 6 curve for a/w ^0.5* 

(ii) J values obtained by all the three models are in 

0 

good agreement with experimental result. The model in 
which all plies fail simultaneously predicts a higher 
value of J than experimental value. Models in which plies 
fail successively predict lower values of • The model 
in which no load relaxation in failed ply is assumed. ,fits 
the experimental results better than the other two models. 



(iii) Difference "between J values for the two models with 
plies failing successively is about 10$, So that the effect 
of load relaxation in failed plies on J is marginal- 
This leads to the conclusion that extensional stiffness 
matrix of tho laminate controls the fracture toughness and 
stacking sequence may not influence it significantly, 

4. 1 SUGGESTIONS FOR FUTURE WORK 

(i) Experimental investigation for DEN and centrally 
located crack specimens should he conducted to find out 

J -integral for those specimens. 

(ii) A three dimensional analysis of laminates may 
be carried out to understand tho effect of interlaminar 
stresses and stacking sequence on J of the laminate. 

(iii) A failure criterion to predict ply failure in a 
cracked specimen should be established, 

(iv) J -integral should be obtained for laminates 
with off axis plies only, c,g. , angle ply laminates (+ e) 
and models presented in tho present study should be 
compared. 
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